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ABSTRACT 

Context. Fossil systems are galaxy associations dominated by a relatively isolated, bright elliptical galaxy, surrounded by a group of 
smaller galaxies lacking L* objects. Because of this extreme environment, fossil groups (FGs) are ideal laboratories to study the mass 
assembly of brightest group galaxies (BGGs). 

Aims. We analyzed the near-infrared photometric and structural properties of a sample of 20 BGGs present in FGs in order to better 
understand their formation mechanisms. This represent the largest sample studied to date. 

Methods. /f ,-band deep images were used to study the structural properties of our sample galaxies. Their surface-brightness distri- 
bution was fitted to a Sersic profile using the GASP2D algorithm. Then, the standard scaling relations were derived for the first time 
for these galaxies and compared with those of normal ellipticals and brightest cluster galaxies in non-fossil systems. 
Results. The BGGs presented in this study represent a subset of the most massive galaxies in the Universe. We found that their 
ellipticity profiles are continuously increasing with the galactocentric radius. Our fossil BCGs follow closely the fundamental plane 
described by normal ellipticals. However, they depart from both the logiTo vs. logLj^^ and log^e vs. logLjf, relations described by 
intermediate mass ellipticals. This occurs in the sense that our BGGs have larger effective radii and smaller velocity dispersions than 
those predicted by these relations. We also found that more elliptical galaxies systematically deviate from the previous relations while 
more rounder object do not. No similar correlation was found with the Sersic index. 

Conclusions. The derived scaling relations can be interpreted in terms of the formation scenario of the BGGs. Because our BGGs 
follow the fundamental plane tilt but they have larger effective radii than expected for intermediate mass ellipticals, we suggest that 
they only went through dissipational mergers in a early stage of their evolution and then assembled the bulk of their mass through 
subsequent dry mergers, contrary to previous claims that BGGs in FGs were formed mainly by the merging of gas-rich galaxies. 

Key words, galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation - galaxies: fundamental parameters - 
galaxies: photometry - galaxies: structure 



1. Introduction 

Fossil groups (FGs) are observationally defined as X-ray bright 
systems (Lx > IO'^^/i^q erg s"'), dominated in the optical by a 
giant elliptical galaxy at the center, and vyith an j^-band differ- 
ence in the absolute magnitude of AMn > 2 dJones et al ] l2003h 
between the two brightest galaxies located within half the virial 
radius of the system. These relatively rare systems were first 
identified by Ponmaji et aL ( 1994) and were the targ et of several 
studie s in the last years ( Vikhlinin et al. 1999; I Jones et alj 
l2003t iMendes de Oliveira et ah 2QQ6i .Cvpriano et al.l l2006t 



Khosroshahi et al ] 120061; iMendes de Oliveira et al. 



20091; 



Zibetti et al.ll2009l;lDemocles et al.ll2010tlAguerri et aLlboi Ih . 

Observationally, these systems are massive galaxy associa- 
tions with typical masses of rich groups or poor clusters. They 
follow also the scaling relations of groups or clusters such as X- 
ray luminosity (Lx) vs. gas temperature (Tx), total mass (M) 
vs. Tx, gas entropy vs. Tx , Lx vs. cluster velocity disper- 
sion (cTf,]), optical luminosity (Lopt) vs. cr^i, 0-^ vs. Tx , ga s 
fraction (/gas) vs. Tx (Khosroshahi et al. 2007; .Sun et al.ll2009h . 
Some differences have been found in the optical vs. X -ray lu- 
minosity (Lopt-Lx) relation dKhosroshahi et al] l2007h even if 
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this is still controversial (IVoevodkin et al.l r2010h . Detailed X- 
ray observations of some FGs also indicate that these systems 
were assembled at early epochs in high centrally concentrated 
dark matter (DM) hal os with large mass-to-light ratios (M/L) 
dDemocles et alj|2010l) . Nevertheless, they do not show cooling 
cores as those detected in galaxy clusters, which points toward 
the presence of other heating mechanisms, like active galac- 
tic nuclei (A GN) feedback (Sunet al 2004; Khosroshahi et al, 
5004 . 2006; Mendes de Oliveira et al.ll2009l) . 

From numerical simulations, and due to the interest- 
ing nature of these systems, their formation scenarios have 
been matter of debate in the last years (D'Onghia et al.' 
2005t ISommer-Larsenll2006!; von Benda-Beckmann et al..,2008i; 



Romeo et alJl2008t iDia z-Gimen ez et alJl201 lUCui et al.ll201 ih 



One scenario attributes their properties to their dynamical his- 
tory, implying that fossil systems are formed if a group/cluster 
remains undisturbed for a long period of time. Therefore, if fos- 
sil systems were assembled in an early epoch of the Universe 
throughout a fast and efficient process of merging, they should 
had time enough for their L* galaxies to merge, thus pro- 
ducing the observed lack of intermediate-luminosity galaxies 
and the large magnitude gap between the bright est and the 
second brightest galaxy of the g roup tD'Onghia et alj 120051; 
Ivon Benda-Beckmann et all l2008h . An alternative explanation 
for FG formation is that such systems are failed groups that 
suffer from a la ck of L* galaxies as an accident of birth 
dMulchaev & Zabl udoff 1999). On the other hand, some authors 
suggest that FGs might represent a transitory phase in the clus- 
ter life, with an absence of significant mergers for a long time , 
enough for clus ter relaxation Avon Benda-Beckmann et alj2008t 
ICuietal.l l2011h. 

The merger history underlying the origin of the central galax- 
ies in FGs continue to be discussed in the literature. Their bright- 
est group galaxies (BGGs) are among the most massive galaxies 
known in the Universe. Originally, they were thought to be pos- 
sible results of group evolution as driv en by dynamical friction 
(iPonman et alJll994tljones et alJl2003h . Recent N-body simula- 
tions and semi-analytic calculations have suggested that BGGs 
formed in a short time-scale as consequence of mergers with low 
angular momentum (ISomme r-Larsen 2006), pushing dynamical 
friction to a less likely scenario for the evolution of BGGs. 

Observations show that BGGs have also different observa- 
tional properties than other bright elliptical galaxies. In partic- 
ular, they present disky isophotes in the center and their lumi- 
nosities correlate with the velocity dispersions of the groups. 
These different properties suggest a different formation sce- 
nario for bright elliptica ls in fossil and non-fossil systems 
dKhosroshahi et alj 1200 6*). Whereas bright ellipticals in FGs 
would grow out of gas-rich mergers, giant ellipticals in non- 
fossil systems would suff er more dry mergers. H owever, neither 
recent samples of BGGs ("La Barbera et al.l l2009l) nor numerical 
simulations IP iaz-Gimenez et al. 2008) show these differences 
dLa Barbera et al.l i2009). All previous results on fossil systems 
have the drawback that they were obtained using small samples 
of FGs. This could be the reason of some contradictory findings 
found by different studies. The lack of a large and homogeneous 
statistical study of this kind of systems makes the previous re- 
sults not conclusive. 

In t he framework of th e Fossil Groups Origins (FOGO) 
project dAguerri et al.ll20ll]) . which aims to carry out a system- 
atic, multiwavelength study of a sample of 34 FGs selected from 
the Sloan Digital Sky Survey (SDSS; Santos et al. 2007), near- 
infrared images in the /Tj-band were taken for the central re- 
gions of 17 groups using the Long-slit Intermediate Resolution 



Infrared Spectrograph (LIRIS) at the 4.2-m William Herschel 
Telescope (WHT). This represent the largest sample of BGGs 
studied to date at these wavelengths. Near-infrared observations 
have the advantage of mapping the distribution of the mass- 
carrying evolved stars and diminish the influence of dust. In 
this work, we analyze in detail the photometric properties of 
the BGGs present on these FGs to perform a systematic study 
of their structural properties and shed light on their forma- 
tion/evolution scenarios. 

The remainder of this paper is as follows: the data sam- 
ple, observations, and data reduction are summarized in Sect. |2] 
In Sect. [3] we describe the surface photometry and photometric 
decomposition. The scaling relations are presented in Sect. |4] 
The discussion of the results and conclusions are presented in 
Sects. |5j and|6l respectively. 

Throughout this paper we assume a flat cosmology with 
Q.TSA - 0.3, Qa = 0.7, and a Hubble parameter Hq = 70 km 
s"' Mpc"'. 



2. Data sample, observations, and data reduction 

2.1. LIRIS data 

Near- infrared imaging of a sample of 17 FG systems belong- 
ing to the FOGO project was carried out at the WHT in 
January 2010. We observed in the /T^-band using the LIRIS 
imaging mode. LIRIS is a near-infrared (0.9-2.4 fj.m) instru- 
ment with an optical system based on a classical camera design 
(Manchado et al. 2004). The detector is a Hawaii 1024x1024 
HgCdTe array operating at 70 K with a pixel scale of Of.'lS 
pixel ^ 

We exposed a total of 3000 s per field, developing several cy- 
cles of a 5 points dithering for each target. The exposure time of 
individual frames was 10-12 s. The dithering cycling procedure 
allowed us to estimate and then correct the background contri- 
butions of the sky and its significant oscillations. In addition, the 
dithering and flat frames helped us to clean cosmetic effects of 
the CCDs. After these corrections, we obtained maximum vari- 
ations of about 0.4% in the background across the full image. 

Data reduction was carried out under IRAlQ environment 
and LIRIS specific packages. In order to keep the quality of the 
point spread function (PSF) across the whole field of view after 
stacking individual frames, the geometry effect of the distorted 
field was corrected. In the combined images the PSF is consis- 
tent with that of the single ones. This is true even close to the 
edges of the field. 

Our observations were carried out under photometric condi- 
tions and within a seeing range of O'.'6-l'.'l (FWHM) as mea- 
sured by fitting a Moffat PSF to field stars in our images. 
The photometric calibration was performed by observing United 
Kingd om Infra -Red Telescope (UKIRT) faint near-infrared stan- 
dards dCasali & Hawarden 1992) . We have checked our pho- 
tometric calibration taking the Two Micron All Sky Survey 
(2MASS) photometry as reference. Our zero-point estimations 
were in agreement with those of 2MASS within ±0.02 mag. 

Finally we identified galaxies in our Ks-hand images 
and measured their ma gnitudes with the SExtractor package 
dBertin & Arnoutslll996 ). Objects were identified imposing that 
they cover a certain minimum area and have a number counts 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1: /Tj-band counts per square degree per 0.5 magnitude bin 
for the fields observed with LIRIS. The continuous curve and 
vertical dotted lines represent the logarithmic function fit to the 
histogram, completeness, and limit magnitudes, respectively. 



above a limiting threshold taking the local sky background as a 
reference. The limiting size and flux were 9 pixels and 1 .5 stan- 
dard deviation of the sky counts, respectively. The selected lim- 
iting size corresponds to an apparent size of (X'TS, which is about 
the typical seeing FWHM obtained from our combined images. 
We performed careful visual inspections of the frames in order 
to deal with the best combination of the above parameters that 
removes spurious objects from the catalogues. 

Figure 1 shows the /Tj-band counts for the fields observed 
with LIRIS. The counts follow the count-magnitude relation ex- 
pected for a homogeneous galaxy distribution in a universe with 
Euclidean geometry. The completeness was computed as the 
magnitude where the measured counts are less than 10% with 
respect to the fit, whereas the limit magnitude was assumed to 
be that of the faintest detected objects. For our sample, the com- 
pleteness and limit magnitudes correspond to targets with signal- 
to-noise ratio S /N=5 and 3, and are Kg = 18.7 and 20.3, respec- 
tively. The characteristics of the observed sample are listed in 
Table B 

2.2. Additional data 
2.2.1. UKIDSS data 

Since we were not able to observe the whole sample of 34 FGs of 
the FOGO sample, we carried out a search in the UKIRT Infrared 
Deep Sky Survey (UKIDSS) eighth data release database to 
complete our sample. We found 7 galaxies out of 17 already 
observed with LIRIS. In addition, 5 new groups not previously 
observed with LIRIS were found. 

Since the quality and depth of the UKIDSS data is clearly 
poorer than those of images obtained using LIRIS, we analyzed 
the images of the 7 common BGGs (see Table. [T]) to under- 
stand whether the UKIDSS data could be useful to our scien- 
tific purpose. We perform a photometric decomposition of the 
surface-brightness distribution of the BGGs measured on the 
UKIDSS images as we did for the LIRIS data (see Sect. [TT] 
for details). By comparing the total magnitudes obtained from 
both fits, we found that only for 4 BGGs (see Table. [T]) the re- 
sults are comparable within 0.2 mag. After analyzing the re- 
sults of the fits we concluded that in order to obtain a reli- 



able fit we need a suitable combination of surface-brightness 
dynamical range and Sersic index, i.e., we found that galax- 
ies with Sersic index « > 3 need at least a range of 6 mag 
arcsec"^ to be well fitted within the observational errors; instead, 
for n < 3 we need only a range of 4 mag arcsec"^. Applying 
these criteria to the 5 images of the BGGs without counterpart in 
LIRIS we found that three groups (SDSSJ154855.85H-085044.3, 
SDSSJ16143L10H-264350.3, SDSSJ225630.04-003210.8) sat- 
isfy the conditions and only these were included in the analysis. 

It is worth notice that we are not covering the complete sam- 
ple of 34 FGs comprised by the FOGO project. This subsample 
of 20 groups is missing mainly the fainter end of the FOGO sam- 
ple in both optical magnitude of the BGG and X-ray luminosity 
of the groups. However, this bias does not affect the results ob- 
tained in this article. In Fig. |2]we show the distribution of the 
complete FOGO sample and that studied in this paper as func- 
tion of the redshift. X-ray luminosity, and r-band absolute mag- 
nitude. 

2.2.2. SDSS data 

Since all our BGGs were observed spectroscopically by the 
Sloan Digital Sky Survey (SDSS), we searched for the presence 
of velocity dispersion measurements in the SDSS database. We 
found the stellar velocity dispersi ons for the whole sa mple using 
the SDSS-III eighth data release (lAihara et al.ll2011l) . However, 
since the SDSS spectroscopy was performed using fibers with 3" 
of diameter, the measurement of the velocity dispersion could be 
contaminated by the presence of another objects inside this ra- 
dius. For this reason, we checked visually the galaxies in order 
to take into account for this effect and we did not use the velocity 
dispersions of the BGGs in the SDSSJ084257.55H-362159.2, and 
SDSSJ133559.98-033129.1 groups. Therefore, our final sample 
of velocity dispersion includes 18 galaxies. The velocity dis- 
persions were corr ected for the aperture size effect following 
the prescription bv lJ0rgensen et al.1 (1 19951) and adopting a fixed 
physical aperture of re/8. The measurements after aperture cor- 
rections are given in Table [1] 

3. Surface photometry 

3. 1. Pliotometric decomposition 

The two-dimensional photometric decomposition of the sam- 
ple galaxi es was performed us i ng the GASP2D algorithm de- 
scribed in lMendez-Abreu et al. I (l2008l) . This algorithm has been 
successfully applied to several galaxy samples jMo relli et aL 
'2008; Pizzella et al. 2008; Beifiori et al. 201 U lPalkBonta et al. 
.201 1) . In particular, it has been applied to a large sample 
of brightest cluste r galaxies (BCGs) in non-fossil systems 
(lAscasoetal.ll201lh 

The GASP2D algorithm is based on a minimization. 
Therefore, the choice of the initial trials for free parame- 
ters is crucial to obtain a good fit. To this aim, the ellipse- 
averaged radial profiles of surface brightness, ellipticity, and po- 
sition angle were analyzed by following the prescriptions by 
Mendez-Abreu et al. (2008). Starting from these initial trial pa- 
rameters the different photometric models of the surface bright- 
ness were fitted iteratively by GASP2D to the pixels of the 
galaxy image to derive its photometric structural parameters. 
Each image pixel was weighted according to the variance of 
its total observed photon counts due to the contribution of both 
galaxy and sky, and determined assuming photon noise limita- 
tion and taking the detector read-out noise into account. Seeing 
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Fig. 2: Distribution of redshift (left panel), X-ray luminosity (middle panel), and r-band absolute magnitude (right panel) for the 
FOGO FGs (solid black line) and of FGs studied in this work (red dotted line). 



Table 1: Characteristics of our sample of FGs and BGGs. 



Group 




z 


m,- 




Source 


FWHM 




(xlO'*' erg/s) 




(mag) 


(km/s) 




(arcsec) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SDSSJ015021 .27- 100530.5 


14.40 


0.365 


17.26 


341+44 


L 


0.7 


SDSSJ015241.95+010025.5 


15.10 


0.230 


15.72 


317±16 


L,U 


1.1 


SDSSJ075244. 19+455657.3 


0.562 


0.052 


14.46 


221 ±5 


L 


0.8 


SDSS J080730.75+34004 1 .6 


4.210 


0.208 


16.38 


263+21 


L 


0.7 


SDSSJ084257.55+362159.2 


29.50 


0.282 


16.79 




L 


0.8 


SDSSJ084449.07+425642. 1 


0.211 


0.054 


14.08 


157±4 


L 


0.7 


SDSSJ090303. 18+273929.3 


17.40 


0.489 


18.06 


273±54 


L,U* 


0.7 


SDSSJ094829.04+495506.7 


6.260 


0.409 


18.21 


270±30 


L 


0.6 


SDSSJ104302.57+005418.2 


4.990 


0.125 


15.98 


235±8 


L,U* 


0.9 


SDSSJ105452.03+5521 12.5 


11.70 


0.468 


17.69 


365+61 


L 


0.7 


SDSSJl 1 1439.76+403735. 1 


4.190 


0.202 


17.14 


218+11 


L 


0.6 


SDSSJl 12155.27+104923.2 


3.850 


0.240 


16.97 


251+21 


L,U* 


0.7 


SDSSJl 14128.29+055829.5 


2.190 


0.188 


16.03 


309±28 


L,U 


0.8 


SDSSJl 14647.57+095228.1 


4.930 


0.221 


16.36 


276+18 


L,U* 


0.6 


SDSS J 1 24742.07+4 1 3 1 37.6 


0.625 


0.155 


15.88 


240±13 


L 


0.6 


SDSSJ133559.98-033129.1 


3.680 


0.177 


15.84 




L 


0.9 


SDSSJ154855. 85+085044.3 


0.509 


0.072 


13.50 


326±9 


U 


0.7 


SDSSJ161431. 10+264350.3 


2.370 


0.184 


15.76 


276±15 


U 


0.7 


SDSSJ225630.04-003210.8 


2.180 


0.224 


16.81 


274+20 


u 


0.7 


SDSSJ235815. 10+150543.5 


0.926 


0.178 


16.08 


297+14 


L,U 


0.8 



Notes. Col. (1): Group name; Col. (2): X-ray luminosity from lSantos et al.l ( 120071) : Cols. (3, 4, 5): redshift, r-band model magnitude, and aperture 
corrected velocity dispersion of the BGGs derived from SDSS; Col. (6): source of the photometric observations: L=LIRIS, U=UKIDSS, common 
galaxies with successfully fitted UKIDSS images are marked with an asterisk; Col. (7): FWHM of the Moffat PSF fitted to the stars in the image 
field of view. 



effects were also taken into account by convolving the model 
image with a circular Moffat PSF with a FWHM matching the 
observed one (see Table [T]i. The convolution was performed as 



a product in Fourier domain before the least-squares minimiza- 
tion. 

We tested three different parametric models to fit the 
galaxy surface-brightness distribution. The de Vaucouleurs pro- 
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file (Ide Vaucouleursl[r948h was extensively used to describe the 
light distribution of elliptical galaxies and to look for central 
dominant galaxies (cDs) in clusters by searching light excess 
over th e model de Vaucouleurs profile in the o uter parts of the 
profi le jMatthews et al.ll 19641: ISchombertH"987l) . The Sersic pro- 
file (H ersiq is nowadays believed to be the most suit- 
able model to describe the surface-brightness distribution of el- 
liptic al galaxies in a wide rang e of lum inosities and masses 
(e.g., ICaon et al. 1993; Graham & Guzman 2003; Aguerri et al.l 
l2004t iKormendv et al...2009) . It is parametrized by re, le, and 
n which are respectively the effective radius, surface brightness 
at re and a shape parameter describing the curvature of the ra- 
dial profile. When n = 4 we recover the de Vaucouleurs pro- 
file and for n = 1 we obtain the exponential one. Figure [3] 
shows the best fit of the observed surface-brightness radial pro- 
file to a Sersic model obtained for our sample galaxies. In ad- 
dition to the previous models, we also fitted a two-component 
(SersicH-exponential) model to the data. This allowed us to fit 
the low surfac e-brightness and extended component that BCGs 
could present ('Nelson et al."2002'; 'Gonzalez et al.' '2^031 120051: 
[Seigar et al. 2007; Liu et al. 2008; Ascaso et al. 201 li We found 
that 7 galaxies are described by this two-component model, how- 
ever, since all galaxies can be successfully modelled by a Sersic 
profile and in order to work with a self-comparable set of struc- 
tural parameters, we decided to use only the Sersic structural pa- 
rameters in this paper In addition, we found that a high fraction 
(14 out of 20) of our sample galaxies are compatible within the 
errors bars with follo wing a de Vaucou leurs profile. Following 
the criteria given by ISchomber3 (Il987h to identify cDs based 
on looking for an upward break in the surface-brightness profile 
with respect to the typical de Vaucouleurs profile, we identified 
SDSSJ124742.07H-413137.6 as the only one candidate to be a 
cD galaxy. However, this classification should be interpreted as 
tentative, since the surface-brightness limits reached in our pho- 
tometr y are not directly comparable with those of ISchombertI 
([T987h . 

The parameters derived for the structural components of the 
sample galaxies using a Sersic fit are collected in Table |2l The 
values of this table were corrected for galaxy inclination, cos- 
mological dimming, K-correction, and evolutionary effects. K- 
correction an d evolutionary effe ct were taken into account using 
the recipes bv iPoggiantil d 1 997h . All the quantities plotted in the 
figures of this paper show such corrected values, unless other- 
wise stated. 



3.2. Internal errors in the structural parameters 

The errors given in Table |2] for every structural parameter were 
obtained through a series of Monte Carlo simulations. Due 
to the fact that formal errors obtained from the minimiza- 
tion method are usually no t representative of the real errors 
dMendez-Abreu et all lloosh . we carried out extensive simula- 
tions on artificial galaxies in order to give a reliable estimation of 
these errors. Monte Carlo simulations have also the advantage of 
getting rid of systematic biases in our measurements, as for ex- 
amp le those introduced by the effects of seeing in our images 
(see lTruiillo et alJl200lh . 

A set of 2000 images of galaxies modeled with a Sersic 
profile was generated. The structural parameters of the artificial 
galaxies were randomly chosen among the following ranges: 



In order to mimic the same instrumental setup, we added a 
background level and photon noise to these artificial images sim- 
ilar to the observed images. They were also convolved simulat- 
ing the seeing that we have in our observations. Finally, these 
simulated galaxies were fitted with the same conditions as the 
real ones. Then, simulated galaxies were used to determine the 
errors of the fitted structural parameters. To assign to every sin- 
gle galaxy the corresponding error for every structural param- 
eter, we divided our catalogue of artificial galaxies in bins of 
0.5 mag, we assumed that the errors were normally distributed, 
with mean and standard deviation corresponding to the system- 
atic and typical error, respectively. Then, we placed our observed 
galaxy in its magnitude bin and assigned to every parameter the 
corresponding error. 

3.3. Radial profiles 

The surface photometry was measured by fitting ellipses to the 
isophotes of the galaxies using the IRAF task ELLIPSE. This 
routine fits a large number of free parame ters to every i s ophot e 
using the iterative method described by iJedrzeiewskil (Il987h . 
Bad pixels, foreground stars, and other spurious sources were 
masked before performing the fit. Figure |4] shows the elliptic- 
ity (e), centroid position ix^),y()), fourth order cosine Fourier 
moment (04), and position angle (PA) radial profiles of the 
isophotes. It is worth noticing that the surface-brightness pro- 
files of 16 galaxies extend at least to 2re at a level of fi^s ~ 21 
mag arcsec"^ corresponding to 1 cr of the sky counts, and 9 reach 
at least 3re, denoting thus the depth of our images. 

No galaxies with ellipticities larger than 0.7 were found, 
in agreement with typical ellipticals (Fig. |4]\). We found that 
in almost all cases the ellipticity profile increases with radius 
up to large distances. We checked that this effect was neither 
caused by a bad sky subtraction nor by light contamination from 
a close object. This behaviour cannot be explained by assuming 
that these galaxies are oblate (or prolate) spheroids and therefore 
suggests the presence of another component or a structural dis- 
tortion in the outer parts (Porter et al. 1991: Men dez-Abreu et alJ 
2010) . There is no correlation with the Sersic shape parameter, 
whereas we found a weak dependence of the ellipticity with the 
galaxy magnitude. Less luminous galaxies are rounder 

The variations of the galaxy centroid do not show any cor- 
relation with the magnitude (Fig. |4)3). Nevertheless, these vari- 
ations seem to be larger for radius greater than re at least for 
some galaxies, possibly indicating the presence of distortions in 
the outer parts of these galaxies. 

The fourth order cosine Fourier coefficient is related with the 
isophotes shape. Negative values indicate boxy isophotes. On the 
contrary, positive values are related with disky isophotes. For our 
sample galaxies these profiles are too noisy and this prevented 
us to assign a single value, or even define a trend (Fig. |4f ). 
Therefore, no obvious correlation was found between this coeffi- 
cient and either the Sersic shape parameter or galaxy magnitude. 

The position angle profiles are shown in Fig. |4lD. We did 
not find any correlation between these profiles and other galaxy 
parameters. In addition, the profiles are very flat with variations 
I PA - PA(re) |< 10°. 



4. Scaling relations 



1 < re < 20 kpc; 
0.5 < n < 6; 



0.5 < ^ < 0.9 

1Q<K,< 14 mag 



We studied the structural scaling relations for our sample of 
BGGs in FGs by analyzing the results from the photometric 
(1) decomposition presented in Sect. 13. ll We also considered the 
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Fig. 3: Observed (black points) and model (green solid line) surface-brightness radial profiles of the sample galaxies. The observed 
radial profiles were obtained by fitting ellipses to the isophotes of the /Tj-band images. The model profiles follow a Sersic law. These 
profiles were not corrected for cosmological dimming, K- or evolutionary corrections. The short vertical lines mark the effective 
radius. Semi-major axes of fitted ellipses are given in arcsec (bottom) and kpc (top). The inset in the upper panels shows the surface 
brightness in logarithmic scale to enhance the quality of the fit in the inner regions. The residuals of the fit are represented with a 
grey region in the lower panels. 



sample of early-type galaxie s by Pahre et al.l (Il998l hereafter 
P98), described in detail in iPahre ( 19991) . as a control sam- 
ple. This represents a statistically significant sample of galax- 
ies in clusters and groups observed homogeneously in the K, 
band. We selected the P98 BCGs to compare the same type 
of objects. The P98 BCGs showed in the figures are the fol- 



lowings: NGC 4874 and NGC 4889 (Coma), NGC 545 and 
NGC 547 (Abell 194), NGC 6166 (Abell 2199), NGC 7720 
(Abell 2634), NGC 4696 (Centaurus), NGC 1316 (Fornax), 
NGC 3309 and NGC 3311 (Hydra), NGC 1272 and NGC 1275 
(Perseus), NGC 4486 (Virgo), and NGC 1407 (Eridanus). 
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Fig. 3: Continued. 
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4. 1. Consistency with P98 data 

The P98 data represents a large sample of early-type galaxies 
observed in the near-infrared with a photometric depth similar 
to that of our images. Therefore, it is ideal to compare with our 
sample of BGGs in fossil systems. However, their structural pa- 
rameters {Lk^, re, and (yUg)) were derived from circular aperture 
photometry which differs from our approach of performing two- 
dimensional photometric decompositions. 

In order to check for possible biases between the two datasets 
that could contaminate our results, we calculate the structural 



parameters for our sample data using also circular aperture pho- 
tometry onto our images. The actual computation was done us- 
ing the IRAF package ELLIPSE, as done by P98. Figure |5] 
shows the results from this comparison. A clear bias appears in 
all the photometric parameters under study. The luminosity (La,), 
effective radius (rg), and mean effective surface-brightness ({fie}) 
derived by the photometric decomposition are 0.11 dex, 0.13 
dex , and 0.44 mag arcsec"^ larger than those derived by aper- 
ture photometry, respectively. These differences are not unex- 
pected and have bee n poi nted out in other studies. F or instance, 
lAguerri et all (l2005h and iTruiillo & Aguerril (|2004 found that 
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Fig. 3: Continued. 



the differences in magnitudes between the two methods is never 
larger that 0.5 mag, in agreement with our result. The main 
reason for these differences is that models used in the photo- 
metric decomposition approach are extrapolated to the infinite, 
therefore given larger magnitudes, effective radius and effective 
surface-brightness. 

Figure |5] also shows the comparison between the central ve- 
locity dispersion of the galaxies using P98 and this work. As 
explained in Sect. 12.2.21 we corrected our velocity dispersions 
using a physical aperture of re/8 whereas P98 used a physical 
scale of 1.53 kpc. These different aperture corrections do not 



introduce a bias between the two samples and it only could in- 
crease the scatter in the relations. 



These biases must be kept in mind in the following but they 
are not enough to change the results presented in this work. In 
fact, most of the differences are within the error bars of the mea- 
surements. To avoid confusions, we have corrected for these bi- 
ases in all the figures where we compare our measurements with 
P98 data. 
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Fig. 3: Continued. 



Table 2: Structural parameters in the /Tj-band derived by fitting a Sersic model to our sample of BGGs 



Group 








n 


q 


PA 






(mag/arc sec^) 


(arcsec) 


(kpc) 






(°) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



SDSSJ015021.27-100530.5 19.16±0.24 4.14±0.73 21.0±3.7 3.68±0.40 0.72+0.02 144.4±4.4 14.08+0.28 

SDSSJ015241.95+010025.5 19.50±0.12 6.77±0.51 24.9±1.9 2.39+0.12 0.52+0.01 18.5+0.2 13.12+0.10 

SDSSJ075244. 19+455657.3 17.27+0.07 3.43+0.14 3.5 +0.1 4.56+0.10 0.95+0.01 125.4+1.4 11.36+0.09 

SDSSJ080730.75+340041.6 20.23+0.39 8.67+2.75 29.5+9.3 3.57+0.48 0.80+0.02 73.4+1.5 13.03+0.32 

SDSSJ084257.55+362159.2 18.39+0.42 3.51+1.09 15.0+4.6 3.59+0.73 0.82+0.02 89.7+1.2 13.41+0.33 

SDSSJ084449.07+425642.1 18.19+0.15 6.05+0.54 6.4 +0.6 3.84+0.21 0.70+0.01 20.9+0.1 11.15+0.08 

SDSSJ090303. 18+273929.3 19.14+0.38 3.56+1.29 21.5+7.8 2.02+0.43 0.73+0.05 119.9+5.3 15.08+0.47 

SDSSJ094829.04+495506.7 21.00+0.42 8.81+3.09 48.0+16.8 5.33+0.39 0.45+0.01 53.0+1.1 14.22+0.34 

SDSSJ104302.57+005418.2 18.36+0.34 3.11+0.74 7.0+1.7 2.88+0.48 0.76+0.02 151.0+3.0 13.19+0.25 

SDSSJ105452.03+552112.5 18.20+0.25 2.72+0.62 16.0+3.7 1.75+0.22 0.78+0.02 101.3+3.3 14.72+0.35 

SDSSJl 11439.76+403735.1 18.89+0.21 2.72+0.41 9.0+1.4 4.41+0.39 0.93+0.03 82.3+2.5 14.08+0.27 

SDSSJl 12155.27+104923.2 18.61+0.34 2.92+0.79 11.1+3.0 2.71+0.41 0.80+0.03 100.9+3.0 14.03+0.34 

SDSSJl 14128.29+055829.5 19.20+0.37 5.10+1.44 16.0+4.5 2.24+0.44 0.76+0.02 100.7+2.1 13.32+0.31 

SDSSJl 14647.57+095228.1 18.95+0.15 4.58+0.49 16.3+1.7 2.04+0.17 0.54+0.01 55.0+1.1 13.47+0.19 

SDSSJ124742.07+413137.6 19.97+0.46 9.89+3.81 26.6+10.2 3.03+0.63 0.58+0.02 149.5+3.2 12.37+0.36 

SDSSJ133559.98-033 129.1 18.16+0.25 3.53+0.55 10.6+1.6 1.99+0.30 0.80+0.02 131.1+2.6 13.10+0.21 

SDSSJ154855. 85+085044.3 18.39+0.03 14.80+0.30 20.3+0.4 2.82+0.03 0.63+0.01 59.8+0.3 10.65+0.06 

SDSSJ161431. 10+264350.3 18.39+0.08 6.37+0.45 19.7+1.4 2.76+0.14 0.78+0.02 20.5+0.4 12.93+0.18 

SDSSJ225630.04-003210.8 19.29+0.08 7.75+0.54 27.9+2.0 3.87+0.19 0.68+0.01 128.5+2.6 13.37+0.18 

SDSSJ235815. 10+150543.5 18.91+0.33 4.91+1.18 14.8+3.6 2.21+0.39 0.69+0.01 97.2+1.9 13.08+0.27 



Notes. Col.(l): Group name; Col. (2): effective surface brightness; Col. (3): effective radius (arcsec); Col. (4): effective radius (kpc); Col. (5): 
Sersic index; Col. (6): minor to major axis ratio; Col (7): position angle; Col. (8): ^j-band total magnitude derived from the Sersic fit. Effective 
surface brightness and magnitudes were corrected for galaxy inclination, cosmological dimming, K-correction and evolutionary effects. 



4.2. Structural parameters vs. luminosity 

Figure |6] shows the relation between the absolute rest-frame 
magnitude in Ks-hand (or equivalently also the /Tj-band lumi- 
nosity in solar units) obtained from the Sersic fit to the surface- 
brightness distribution of the sample galaxies and the mean ef- 
fective surface brightness, or Sersic shape parameter. 

Notice that some of our sample galaxies are among the 
brightest galaxies located in cluster/group environments in the 
Universe. They are brighter than the BCGs present in nearby 
clusters such as Virgo or Coma, which are represented by trian- 
gle symbols. However, they are neither the most concentrated, 
nor posses the highest central surface brightness. In fact, the 
mean effective surface brightness correlates well with the galaxy 



magnitude and follows the same trend of normal ellipticals and 
BCGs. 



The P98 sample does not provide information about the 
shape of the surface-brightness profile, which in our case is rep- 
resented by the Sersic index, and therefore we can not make a 
direct comparison. However, it is worth noticing the wide range 
of Sersic index values obtained for our BGGs, ranging from 1.7 
to 5.3 with a mean value of 3.0. This indicates that the shape of 
the surface-brightness profile can be very different from a typical 
de Vaucouleurs profile. 
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Fig. 4: (A) Ellipticity, (B) centroid (r - yjix- x^p- + {y — yo)^), (C) fourth order Fourier cosine moment, and (D) position angle 
radial profiles of our sample galaxies. All radial profiles are normalized to the effective radius in the abscissa axis. Blue, green and 
red solid lines refers to BGGs with Mk, > -26, -26 < Mk, < -27, and Mk, < -27, respectively. Both the centroid and position 
angle profile were reported to their values at rg. The radial profiles are plotted for radii larger than the corresponding PSF FWHM. 
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Fig. 6: Correlations between (yUe) (upper panel), n (bottom panel) 
and the galaxy absolute magnitude in the /Tj-band. Black large 
dots indicate our sample galaxies. Blue dots represent the P98 
early-type galaxies. The BCGs in the P98 sample are marked by 
blue open triangles. 



4.3. Fundamental plane 

It is well known that the structural parameters of clas- 
sical bulges and elliptical galaxies follow a tight relation 



called the Fundamental Plane (FP; 'Diorgovsk i & Davislll987l : 
[Dressier et al. 1987) among the effective radius (log re), mean 
effective surface brightness ((/ie)) and central stellar velocity dis- 
persion (log ctq). The small intrinsic dispersion of the FP rela- 
tion is a consequence of the virial theorem and demonstrate that 
these systems are in dynamical equilibrium. Possible deviations 
from the FP are related to a variation of the mass-luminosity 
ratio with the luminosity and/or a non-homology in the struc- 
ture of these galaxies throughout the whole range of luminosities 
(Truiillo et al. 2004). The FP has been widely used as a powerful 
tool in measuring galacti c distances and also in studies of galaxy 
formation and evolution ( Kiaergaard et alJ 19931 : lj0rgensen et al.l 



[T^FTsifgt [Kelson et alJll997l) . 



Figure[7]shows the distribution of our sample galaxies in the 
FP relation defined by the sample of P98. Our BGGs follow a FP 
compatible within the errors with that of normal ellipticals and 
BCGs, therefore implying that in the innermost regions (r < r^) 
they are relaxed systems similar to normal ellipticals. 



4.4. Kormendy relation 

The effective radius and mean surface brightness of elliptical 
galaxies are correlated through the so-called Kormendy relation 
(KR; [Korm endv 1977) which represents a projection of the FP. 

Figure [SI shows the position of our BGGs in the KR with re- 
spect to the control sample of P9 8. The fit to the data points was 
done using the FITEXY routine (IPress et alJll992) implemented 
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Galaxy jf Galaxy § 

Fig. 5: (A) luminosity , (B) effective radius, (C) mean effective surface-brightness, and (D) velocity dispersion differences between 
the aperture photometry and photometric decomposition measurements. Abscissa axes represent a number assigned to each galaxy. 
Blue solid line represent the mean difference taking into account all points. Error bars represent the errors Usted in Tables [T] and |2] 




1.53 log(ao) + 0.314 <fi^> -8.24 



Fig. 7: Distribution of our BGGs (large black points) in the edge-on view of the FP defined by the P98 early-type galaxies (small blue 
points). The BCGs of the P98 sample are marked by blue open triangles. The solid line represents the FP best fit to P98 galaxies. 
The dotted lines represent the 3cr deviation from the fit. 
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in IDL0, which estimates the parameters of a straight-line fit tak- 
ing into account for data errors on both variables. The best fitted 
KR to the P98 sample is given by 

(;Ue) = 15.13(+0.38)H-3.34(+0.86)logre. (2) 

Our galaxies are systematically above the best-fit relation ob- 
tained for the P98 sample, even if both relations are similar once 
the errors are taken into account. The relation obtained by fitting 
only our galaxies is 

(yUe) = 13.48(+2.71)H-3.78(+1.51)logre. (3) 

Several works studied the KR in galaxy clusters demon- 
strating that the intrins ic dispersion is approximately 0.4 mag 
arcsec'^ in (ne) (Hoessel et alJ Il987t ISandage & Perelmuted 

ll99lHLa Barbera et al 20031^ This high intrinsic dispersion can 

be interpreted in different ways: by the fact that the KR does 
not consider the third parameter of FP (i.e., velocity dispersion, 
[Ziegler et al. 1999); because of the measurement errors as well 
as the systematic errors due to the photometric calibration and 
also due to the corrections introduced for different biases (e.g., 
zero point and color transformation, K-correction and redden- 
ing); or because the position of the galaxies in the KR depends 
on their magnitude (Nigoche-Netro et al. 2008). 

In the case of our BGGs, we suggest that since our galaxies 
represent the brightest galaxies located in cluster/groups, the dif- 
ferent position of the BGGs in the KR is produced only by their 
higher luminosity. In fact, galaxies with increasing magnitudes 
are placed parallel to the KR but with higher values of the mean 
surface brightness for a given effective radius, as it is the case 
of our sample galaxies. In addition, the coefficients of our fit are 
poorly constrained due to our small number statistics, and both 
samples are consistent within the errors. 

4.5. Faber-Jackson relation 

The Faber-Jackson relation (FJ; iFaber & Jacksoniri976 h repre- 
sents another physically significant projection of the FP. It cor- 
relates the galaxy magnitude with the stellar velocity dispersion 
and for this reason it is observationally more expensive to obtain 
than the KR since it requires for spectroscopic data. For the sake 
of clarity, in Fig.|9]we plot the total luminosity in the A'5-band in- 
stead of the magnitude. The Ks-hand luminosities of a given stel- 
lar population are less affected by metallicity effects than those 
in optical passbands. In addition, they are a good proxy of the 
total stellar mass since the typical M/L ~ 1 for an old stellar 
population (e.g., Bruzual & Charlot.2003.) . 

Figure |9] represents the FJ r elation for the P98 sarn ple and 
our sample of BGGs. Recently, iBernardi et alJ (1201 lallbl) have 
demonstrated that the scaling relations of early-type galaxies 
present a change in their slo pe at a characterist ic stellar mass 
(M. ~ 2 X 10" Mq, see also lTortora et alJl2009h . Interestingly, 
this curvature of the scaling relations is not present when the stel- 
lar mass of the galaxies is replaced by the central velocity dis- 
persion. In particular, they find that galaxies with masses higher 
than 2 xlO" M0 have smaller velocity dispersions and larger 
effective radius for a given mass than their counterparts in the 
mass range 3xlO'"< M./M© < 2x10". 

Our sample galaxies represent a subset of the most mas- 
sive galaxies of the U niverse. Some of them probably fall in 
the Bernardi's sample dBernardi et al.ll201 lalbl) since they have 

^ Interactive Data Language is distributed by ITT Visual Information 
Solutions. 



SDSS information. Unfortunately, we cannot compare directly 
with their work since they have optical photometry whereas we 
have near-infrared photo metry. In Fig. [ 9]we pl ot the FJ fit to the 
range of masses where Bernardi et alJ (1201 Ibl) consider that the 
relation is linear. Most of our BGGs follows the expected curva- 
ture in the lo g ctq vs. log relation , confirming and extending 
the results bv lBernardi et al.l (1201 Ibl) . 

4.6. Effective radius vs. luminosity 

Figure [To] shows the relation between the /Tj-band luminosity in 
solar units obtained from the Sersic fit to the surface brightness 
and the effective radius of the galaxies. 

As done in the pr evious section, we tested the results of 
IBernardi et al.l (1201 Ibl) with our sample of massive BGGs. We 
fitted a straight line in logarithmic scale to the data in the mass 
range 3x10'" < M,/Mq <2x10" and extrapolated such a fit to 
larger masses. From our data, it is clear that most of our BGGs 
follow the expected curvature in the logrg vs. logLj^^ relation, 
i.e., galaxies with masses higher than 2 xlO"M0 have larger 
effective radius for a given mass than their counterparts in the 
mas s range between 3x10'° < M,/Mq < 2x10", in agreement 
with IBernardi etaP (1201 Ibl) . 



5. Discussion 

The BGGs of the fossil systems presented in this study repre- 
sent a subset of the brightest and most massive galaxies in the 
Universe. This can be seen by looking at both their /Tj-band lu- 
minosities and central stellar velocity dispersions. However, they 
closely follow the FP scaling relation defined by normal ellipti- 
cal galaxies and other BCGs, suggesting that they are dynami- 
cally relaxed systems at least in their central regions. 

FGs are believed to be created fast and in an early epoch of 
the U niverse by the efficien t merger of L* galaxies ins ide DM 
halos (iD'Onghia et al.ll2005l: iDfaz-Gimenez et al.ll20Tll) . In this 
scenario, the BGGs observed at low redshift had time to reach 
a virialized status although some remnant of their merger phase 
could still be present. The different relaxation time scales be- 
tween the inner and outer regions imply that the outskirts of 
galaxies may have fossil information about their past formation. 
The continuously increasing ellipticity profile observed in al- 
most al l galaxies of our sam ple could be a hint of such a merger 
history. [Porter et al] (1 19911) reported on this behaviour of the el- 
lipticity profiles in a sample of BCGs, they found that 50% 
of the BCGs in rich clusters present these characteristics and 
associated them with the observed absence of r otation in lumi- 
nous ellipticals (Davies et al. 19831 ICarter et al.. 1985). This re- 
sult might be consistent with our findings if the increasing ellip- 
ticity profiles are corr elated with galaxy luminosity. In addition, 
Bovlan-Kolch in et~an (f2006) demonstrated by using numerical 
simulations that these ellipticity trends can be explained in terms 
of a high fraction of radial collisions, indicating a preferential di- 
rection in the orbits of the mergers. 

Assuming that FGs were formed in an early epoch 
(Dariush et al. 2007; Diaz-Gimenez et al. 2011) and then pas- 
sively evolved for a long time in order to reach a relaxed state, 
it is customary to consider that dissipational mergers with a 
high fraction of gas should be predo minant in the formation 
of these systems ( Dfaz-Gimenez et al. 2008). Recent works by 
Robertson et al. (2006) and Hopkins et al. (2008) have demon- 
strated through simulations and observations, respectively, how 
the observed tilt of the FP can be explained in terms of the 
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Fig. 8: The KR for our BGGs (large black points) and P98 early-type galaxies (small blue points). The BCGs of the P98 sample are 
marked by blue open triangles. The sohd line represents the KR best fit to our sample and the dotted line represents the best KR fit 
to P98 sample. 




Fig. 9: Distribution of our BGGs (red stars and large black points) and P98 early-type galaxies (small blue points) in the log ctq vs. 
logLt^ plane. The BCGs in the P98 sample are marked by blue open triangles. The solid line represents the best fit to the galaxies 
in the luminosity range 3x10'" < LkJLq < 2x10" as done by Bernardi et al. (2011b). Red stars and large black points represent 
BGGs with ellipticities e > 0.3 and e < 0.3, respectively. The bottom panel represents the residuals from the best fit. 
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Fig. 10: Distribution of our BGGs (red stars and large black points) and P98 early-type galaxies (small blue points) in the log vs. 
log Lfc, plane. The BCGs in the P98 sample are marked by blue open triangles. The sol id line represents the best fit to the galaxies 
in the luminosity range 3x10'" < LkJLq < 2x10" as done bv lBernardi et alj (1201 Ibh . Red stars and large black points represent 
BGGs with ellipticities e > 0.3 and e < 0.3, respectively. The bottom panel represents the residuals from the best fit. 



amount of dissipation involved in the formation of elliptical 
galaxies. They conclude that galaxies formed through dissipa- 
tionless mergers perfectly follow the FP defined by the virial 
theorem (re oc crQ(yUe)" '*) whereas ellipticals formed through dis- 
sipational mergers will present a FP tilt depending on the frac- 
tion of gas involved in the merger Our galaxies, as normal el- 
lipticals, are well described by a tilted FP, thus showing some 
indications that most of them could be the remnants of mergers 
with a high gas fraction. In addition, numerical simulations show 
that gas-rich mergers at an early stage of the formation of BGGs 
are need ed to account for the observed properties of elliptical 
galaxies jHernquist et alJll993l:lNaab & Ostrike3l2009h . 

However, numerical simulations of gas-rich mergers also 
predict that the remnants will be round and compact galax- 
ies (i.e., with small effective radius) where the stellar veloc- 
ity dispers i ons s h ould increase (e.g., iKhochfar & Silkl 120061: 
iNaab et alj 120071: iHopkins etal] 12008 '). Recent observational 
works on early type galaxies bv iBern ardi et al. (2011ab) have 
shown a change in the slope of different scaling relations that 
involve the stellar mass. This trend is interpreted by the authors 
in terms of the kind of mergers that have been involved in the 
formation of these systems. Major dissipationless mergers be- 
tween galaxies are mechanisms expected to increase the final 
size of the resulting galaxy but not its central velocity disper- 
sion; however, if minor dry mergers are the predominant mech- 
anism, they a re expected to change both the s ize and veloc- 
ity d ispersion (IBernardi et al J 120 1 la*). Therefore. IBernardi et alJ 
(1201 lb) claimed that the most feasible scenario is that galax- 
ies with M» > 2x10"Mq formed mainly by major dissipa- 
tionless mergers. It is clear that most of our BGGs follow the 



expected curvature in both the log ctq vs. log Lk^ (Fig. ^ and 
log Te vs. log Lk, (Fig. [Tol l relations, conf irming the observa- 
tional trend found by B ernardi et alj (1201 lb ). This interpretation 
contradicts the previous one hinting toward a formation scenario 
where most of our BGGs formed through major gas-rich merger 
( Khosroshahi et al.i,2006) . 

Another commonly used proxy to understand the forma- 
tion mechanisms of early-type galaxies is the Sersic index. 
Numerical simulations predict Sersic indexes for the remnants of 
dry rnergers in the rang e 3 < « < 8 (iGonzalez- Garcfa & Balc ellsl 
I2OO5I: iNaaTfe Truiilloll200"6l) . whereas in dissipational rnergers 
they wo uld be in the range 2 < n < 4 dHopkins et al.l 120091 : 
Kormend v et alJl2009l) . We have tentatively divided our sample 
in galaxies with n < 3 and n >3 in order to understand whether 
the shape of the surface-brightness profiles could be used to dis- 
tinguish among the different formation scenarios. We did not 
find any relation between the departure from the intermediate- 
mass scaling relations and Sersic index. Similarly, several works 
based on numerical simulations found that the shape of the 
isophotes can be also used to trace t he merger history of ellip- 
tical galaxies (iNaab & Burkertll2003HKhochfar & Burkertll2005l: 



lNaabetal.l l20b6V These studies found out that disky isophotes 
are the result of gas-rich mergers whereas boxy isophotes are 
preferentially formed in dissipationles s collisions. Supported by 
these results. IKhosroshahi et alJ(l2006l) claim that BGGs in fossil 
systems were formed through dissipational mergers since they 
show non-boxy isophotes. We measured the isophote shapes by 
means of the 04 coefficient, however, their noisy profiles (Fig. HI 
prevented us to extract conclusions about the merger history of 
our BGGs based on the isophotes shape. 
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Interestingly, we found that the ellipticity turned out to be a 
good indicator of the kind of merger involved in the formation of 
BGGs. Numerical simulations show that the remnants of dissipa- 
tional mergers are rounder on av erage when compared to rem- 
nants of dissipationless mergers jCox et al.l |2004). In fact, the 
ellipticity distribution of dissipational merger remnants peaks at 
about 0.2 which is in sharp contrast to the dissipationless rem- 
nants whose distribution of ellipticities peaks at about 0.4. We 
used this information to separate our BGGs in galaxies with 
e < 0.3 and e > 0.3. We measured the ellipticities at re to mimic 
simulations. The results are shown in Figs. |9] and (TO] It is clear 
that BGGs with larger ellipticities systematically deviate from 
the intermediate-mass scaling relation, whereas rounder galaxies 
are also distributed near to this relation. Therefore, we interpret 
this result as a further hint that dissipationless mergers (e > 0.3) 
are the main processes driven the evolution of BGGs. 

We conclude that BGGs in fossil systems suffered major dis- 
sipational mergers in an early epoch of their formation, but the 
bulk of their mass is assembled latter in subsequent dissipation- 
less mergers that increase the BGGs size. 

6. Conclusions 

We analyzed the near-infrared ^T^-band photometric properties 
of a sample of 20 BGGs observed within the framework of the 
FOGO project ( Aguerri et al. 2011). So far, this represent the 
largest sample of BGGs in fossil systems studied in the near- 
infrared. Here we summarize our results: 

- The structural parameters of the sample g alaxies were de- 
rived using the GASP2D algorithm (Me ndez-Abreu et"al] 
[2008) by fitting a Sersic model to the their surface-brightness 
distribution. This model provides a good description for all 
the galaxies and allows a straightforward comparison with 
the results available in literature. 

- Only one sample BGG is a candidate to be a cD galaxy since 
it shows an upward break in the surface-brightness profile 
with respect to a de Vaucouleurs profile. 

- The ellipticity profile of most of our galaxies is an increasing 
function of the radius. This can be related to the fact that our 
galaxies represent the very bright end of the luminosity func- 
tion. The trend in the ellipticity profile cannot be explained 
in terms of a oblate/prolate spheroid. Together with the varia- 
tion observed in the galaxy centroid and position angle radial 
profiles, it might indicate the presence of preferential direc- 
tion in the orbits of the mergers. 

- We built the /Tj-band scaling relations (FP, KR, and FJ) for 
our sample of BGGs and we compared them with those of 
normal elliptical and BCGs. We did not find any significant 
difference between the FP defined by BGGs and those of 
ellipticals or BCGs. This suggests that the central regions of 
BGGs are relaxed with structure similar to those of normal 
ellipticals. However, we found a change of the slope in the FJ 
and in the log rg vs. log Lk^ relation for our massive galaxies. 

These observational results can be interpreted in terms of the 
formation scenario of the BGGs. Our BGGs sample follows the 
FP of ellipticals and therefore are well reproduced by a tilted 
FP. Therefore, this supports that most of our BGGs probably 
suffered dissipational mergers during their formation. However, 
the curvature found in the scaling relation involving the mass of 
the galaxies can be interpreted as an indication that the stellar 
mass of these systems grew mainly by dissipationless mergers. 
Therefore, we suggest that BGGs in fossil system suffered ma- 
jor dissipational mergers in an early epoch of their formation. 



but the bulk of their mass is assembled latter in subsequent dis- 
sipationless mergers that increase the BGGs size. 
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